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1. Introduction 

Matrix quantum mechanics represents an important analytical tool for the investi- 
gation of one dimensional bosonic string theory (see references in review [1]). Matrix 
formulation of this string theory is based on the correspondence between the feynman 
diagrammatic expansion of the matrix functional integral and the sum over all possible 
discretized (say, triangulated) curved two-dimensional manifolds embedded into one di- 
mensional target space (time). 

It was claimed in ref. [2] that the singlet sector of the quantum mechanical matrix 
model (trivial representation of U(N)) compactified on a time circle of the length {3 rep- 
resents the vortex free sector of one dimensional bosonic string theory. This claim was 
supported by a remarkable self-duality with respect to the change (3 — > af//3 expected 
from the continuum string field theory. 

Non-singlet states are of special interest in the ID matrix theory. It was conjectured 
in [2-3] that they correspond to the vortex-anti-vortex excitations (winding modes) of 
compactified ID bosonic string with the target space coordinate living on a circle of a 
finite length (3. In the statistical mechanics language, they describe the classical XY- 
model on a random (dynamical) 2D manifold. 

In ref. [3] the simplest non-singlet case, (the adjoint representation for "angular" 
variables), was investigated. Using the results of the paper [4] where the planar (N — > oo) 
approximation for the adjoint hamiltonian was developed, the authors of ref. [3] have 
demonstrated that there exists a big gap separating ground states for the singlet and the 
adjoint representations: 

E adj -E sing ~ 7^-| log A| (1.1) 

where A = A c — A is the deviation of the cosmological constant from its critical value. 

It was also argued in ref. [3] that this result can be qualitatively reproduced from the 
continuum Liouville theory assuming that one includes a vortex-anti- vortex pair into the 
system. So, the adjoint representation was conjectured to describe a vortex-anti-vortex 
pair in the presence of 2D quantum gravity. Higher representations were associated with 
bigger numbers of vortices. On the base of these conjectures the authors of ref. [3] gave an 
estimate of the position, f3 c , of the Berezinski-Kostelitz-Thouless phase transition arguing 
that (3kt is finite. The phase (3 > (3kt (low vortex density) describes c = 1 system and, 
for (3 < 13 kt, c = (vortices destroy long order correlations). 
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In this paper we shall use a different technique but the result (1.1) will be confirmed 
and even generalized. Namely, we shall calculate the double scaling limit of the parti- 
tion function for the adjoint representation providing the partition function of vortex- 
anti- vortex pair interacting with the 2D quantum gravity to all genera. We shall explicitly 
demonstrate by direct arguments that for the matrix quantum mechanics in the adjoint rep- 
resentation describes the "vortex-anti-vortex" sector of the Xy-model (sometimes called 
as the model of planar rotators) interacting with 2D quantum gravity. As in the XY-mode\ 
on a regular lattice, knowing the partition function in this sector one can determine 0kt 
precisely. Our estimate (3kt = 47r coincides with the one given in ref. [3]. 

We will propose here a mathematical setup, different from [3], for the analysis of 
nonsinglet sectors in the matrix quantum mechanics. Although we can solve a number of 
matrix models exactly we usually need from the solutions only a small piece of information 
about the continuum (thermodynamic) limit of lattice manifolds dominated by huge Feyn- 
man graphs. In this respect it seems to be very desirable to have in hands some formalism 
which has an advantage of the matrix approach to be manifestly invariant (i.e. not using 
the notion of a coordinate system and non-invariant objects like a metric) but avoiding 
the intermediate lattice stage usual for the matrix models and dealing directly with the 
continuum theory. 

In the matrix language one usually diagonalizes the matrices and tries to formulate 
a model in terms of their eigenvalues integrating explicitly, if it appears to be possible, 
over the "angular" degrees of freedom belonging to a group (U(N) in the case of N x N 
Hermite-an matrices) . For many interesting models of 2d gravity with matter fields having 
the central charge c < 1, this integration can be performed exactly. In the resulting 
effective action for the eigenvalues, not the entire potential of the original matrix model 
is important for the continuum (planar or double scaling) limit but only the vicinity of a 
singularity of the effective action. In the most generic case, this is a quadratic top of the 
potential. The behavior of the eigenvalues near the top usually defines the essential, long 
distance phenomena and everything else in the shape of the potential is responsible only 
for microscopic details, i.e. cut-off, area definition etc. This brings to the mind the idea 
that the matrix models can be formulated in a simplified way where only the essential part 
of the matrix potential will be taken into account. One can hope that, in the continuum 
limit, the whole model might appear to be almost (i.e. up to a cut-off) Gaussian and, 
hence, exactly solvable even though the potential is unstable. 
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One of the most transparent examples realizing this mechanism is the matrix quantum 
mechanics [5], which was shown to correspond to the bosonic string field theory in one- 
dimensional target space [6] and was solved later in the double scaling limit [7-10]. This 
theory can be described by matrix quantum mechanics with an upside-down oscillatorial 
potential. The system can be stabilized by putting infinite walls on cut-off distances ±A 
from the top of the potential for every eigenvalue of the matrix field. 

In this case it is natural to expect that interesting physical quantities can be obtained 
by an analytical continuation of corresponding quantities for the standard matrix oscillator 
with the stable quadratic potential, i.e. by the change u — > iu (where u is the frequency 
of the oscillator). 

Life appears to be more complicated: the standard oscillator has more symmetry 
then the upside-down one and it does not know about any cut-off. Usually to establish a 
correspondence between them the direct analytical continuation is not possible and should 
be completed by a guess about the cut-off dependence. One of the recent results of this 
correspondence was given in ref. [11] for the singlet sector Green function where the 
corresponding planar result of ref. [12] was generalized to the double scaling limit. 

In the present paper we give another example of such correspondence. Our goal will 
be to calculate partition functions for some non-singlet states in the ID model in the 
double scaling limit using partition functions for the ordinary matrix oscillator. We argue 
that this correspondence works at least for a sector of the model, given by the adjoint 
representation of U(iV) and (as we show explicitly) describing single vortex-anti-vortex 
configuration. At the end we discuss a possible generalization of this approach to a D + 1 
dimensional matrix model which might give in the future new analytical tools into our 
hands. 

The paper is organized as follows. 

In section 2, we formulate the matrix model with the compactified time and discuss 
its correspondence with ID string theory and the XY-model. 

In section 3, we give a derivation of Hamiltonians in non-singlet sectors in terms of 
the eigenvalues of the matrix field. The double scaling is considered for which the problem 
can be reformulated as the upside-down matrix oscillator. 

In section 4, a U(N) twisted partition function is introduced which allows us to 
expand the oscillatorial partition function into a sum of partition functions of states in 
given representations of U(N). In this section we also show that the adjoint sector of the 
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ID matrix model corresponds to the insertion of one vortex-anti-vortex pair on the string 
world sheet. 

In section 5, the analytical continuation of the adjoint partition function is suggested 
and, in section 6, physical consequences of the obtained results are considered. 

In section 7, a possible generalization of our approach to the case of D + 1 string 
theory is sketched out using the idea of correspondence between Gaussian matrix models 
with the stable and the upside-down quadratic potentials. 

The section 8 is devoted to conclusion. 

Three appendices contain technical details concerning applications of the group theory 
methods to our model. 



2. Matrix quantum mechanics in periodic time as ID compactified string 
theory 

The matrix model under consideration can be defined by the following functional 
integral for the partition function: 

Z N (p,\) = [ V N2 ^(t)exp-Ntr f dt[\^> 2 + V{^)] (2.1) 

J<p{0)=y(l3) Jo 2 

where (fit) is a hermitian-matrix- valued periodic ID field : tpij(t) = Tpj^t); (fij(0) = <fij(P) 
i,j = l,...,N. The potential can be chosen, for example, in the form 

V(<p) = \? - \v\ (2-2) 

which corresponds to triangulations of the string world sheet. 

The standard arguments leading to the string picture interpretation of eq. (2.1) are the 
following [3]. Let us expand logZ^^.X) in the coupling A using diagrammatic Feynman 
rules. The corresponding free energy /tv(/3, A) = log Z N {(3 : A) has the form: 

oo oo „/3 „/J 

f N (P,X) = J2N- 29 J2 XV J2 •••/ dh-.-dtyX 

g=0 V=0 n (V) J° J° 

^9 



IT 5Z e- | * i -*>' +/3mi >' 1 (2.3) 

<ij> mij — — oo 



where g is the genus of a </? 3 -Feynman graph Gg ; V is the number of vertices, n<ij> ^ s 



the product over all links of the graph, < ij >G G g n \ and the periodic propagator in our 



case is 



D(U-tj)= J2 e-^-^ +m ^ (2.4) 

m= — oo 

where t^s are attached to vertices of a graph. The constant A plays the role of a bare 
cosmological constant and 1/N is the (topological) string coupling constant. For sufficiently 
large (3 we expect that only the term with m = in eq. (2.4) will survive and we will come 
to the well investigated case of the ID matrix model on the infinite time interval [1]. In 
this case time coordinates ti,i = 1,2 .. .V, will correspond to the target space coordinates 
of the string world sheet represented by a </? 3 -graph Gg . The corresponding action is of 
the form 

s = ^2 i ti — tj i 



<ij>€G?) 



(2.5) 



which differs from the discretized Polyakov string action 



but is proven to be in the same universality class due to the superconvergency of underlying 
Feynman diagrams. 

Eq. (2.3) describes the XF-model coupled to 2d gravity and we expect the appearance 
of all phenomena related to the vortex configurations of the t-field. The analysis of the 
vortex kinematics in the model was done in ref. [2]. Let us briefly mention it here for 
completeness. Following the classical papers [13], [14], [15] we can view rriij in eq.(2.3) as 
an abelian gauge field defined on links of a </? 3 -graph. The sum along a non-self-intersecting 
loop L on the graph 

(2.7) 



Mi. Yl 1,1 ' 



<ij>EL 

gives the integral charge of the vortices enveloped by the contour L. In the case when L 
is just the boundary of a face of the graph, Ml can be considered as the elementary field 
strength or the vortex number in this face. The duality transformation corresponding to 
the Fourier transform of the original propagators 

oo „ oo 

m=-oo r k i:j = -oo 1 ' V (3 ^JJ 



brings eq. (2.3) to the form where M L are new dynamical variables. If we insert eq. (2.8) 
into eq. (2.3) and integrate over t^s, we obtain the constraint at every vertex i, which is 
just the ordinary condition of the momentum conservation, 



It means that (excluding the zero modes) we have got E — V + 1 = F — 1 + 2g variables 
instead of V — 1 (in virtue of the Euler theorem: F — E + V = 2 — 2g, where F, E and V are 
the numbers of faces, edges and vertices of a graph). The usual choice of these variables is 
to attach a momentum pi to each face and to define remaining 2g variables as momenta 
l a running along 2g independent non-contractable loops on a graph: L a , a = 1, . . . , 2g. 
The resulting dual representation for the partition function (2.3) takes the form 

(n\ 2 oo / j. TO \ —2g / -. I x F +00 +00 

E ? E E - E 

' g=0 \ ' F V 7 MF) P! = -OC p F = -00 

+ OO +OO _ 

y ... y n -2 — 5 (2-io) 

where J2mf) is the sum over all dual graphs (triangulations) having F dual vertices and 
a genus g; Yl<u> ^ s the product over all dual links < I J >G Gg'; e a u = ±1 when a dual 
edge < IJ > crosses an edge belonging to a chosen in advance non-contractable cycle L a: 
and is zero otherwise. The sign has to be chosen with respect to the mutual orientation of 
the link < I J > and the loop L a . 

It can be argued that one can discard the vortices in eq.(2.3) by imposing the "pure 
gauge" conditions on the field (up to a contribution of non-contractable loops) 

rriij = y + rrii- mj (2.11) 

a=l 

where integers l a are defined on non-contractable dual loops L a ; efj is the object dual to 
e a u . Then, eq. (2.10) will look as 
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/oo F oo +oo 2 

■°°J=1 Ii = -oo l 2g = -oo<IJ> 1+ {f ) (PI -PJ + Ea=l l a^j) 
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The sums over pj in eq. (2.10) are substituted by the integrals, since, in the original 
expression (2.3) with the choice (2.11), we can instead of summing over m^'s spread the 
integrations over t^s to the infinite interval. Then the corresponding Fourier transforms 
of propagators will be integrals rather than sums. 

Now the original representation of /jvO#, A), eq. (2.3), with the condition (2.11) looks 
very similar to its dual transform (2.12). The arguments of the propagators and the string 
coupling constant can be matched by the simple duality transformation 



JL 

27 ~J 
1 2tt 



(2.13) 



N (3N 

Two main differences are : the original </? 3 -graphs are substituted by dual ones, and 
the propagators are different. If we believe in universality in the continuum limit , these 
two differences should be insignificant on the macroscopic scale. More important is the 
equivalence of the sets of variables (ti,l a ) and (pi,l a ) in eqs. (2.3) and (2.12). So, one 
can hope that, in the double scaling limit, all answers will be invariant under the duality 
transformation (2.13). It was claimed in ref. [2] that this is really true for the singlet 
sector of the model, which was identified therefore with the vortex-free partition function 
(2.12). We will show in this paper that the inclusion of the non-singlet states destroys the 
self-duality, which is natural to expect in the presence of vortices. 



3. Effective action for eigenvalues 

The first standard step to calculate the partition function (2.1) is the diagonalization 
of the matrix field: 

JV 

^ J (t) = ^0+(t)^(t)O fcj (t) (3.1) 
k=i 

where Q(t) G U(N). 

In terms of new variables we have 

N 

trtp 2 = Y,zf + Z>i " z 3? I A v I' ( 3 - 2 ) 

i=l i^j 

where we introduced "the connection" 

Aij(t) = (Sl+tyij (3.3) 
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The periodicity condition <fiij(t) = <fij(t + (3) implies that eigenvalues, z^, have to be 
periodic in t only up to an arbitrary substitution V: 

N 

Zk{t + (3) = Y J 'Pk j z j {t)V-t 

j=i (3-4) 

which shows that the variables Aij(t) are not independent, contrary to the free boundary 
conditions case, and that they obey the constraint: 

Texp % \ dtAit) = V~ l (3.5) 
Jo 

This makes the integral over A(t) rather complicated unlike the most of solvable mod- 
els of 2d gravity with c < 1. Usually, attempts to solve some matrix models describing 
gravity with matter having the central charge c > 1 are blocked by the difficulties to 
compute the emerging integrals over the "angular" variables. We hope that the investi- 
gation of the model containing at least one nontrivial integration might be useful for the 
understanding of the role of "angular" variables for c > 1. 

The Dyson measure for the new variables z and O looks as 

N 

V<p(t) = n A 2 (z(t))l[dz k (t)d N2 A(t)x 
te[o,/3] fc=i 

x5 u{N) (Texp Adt.V' 1 ) (3.6) 
Jo 

where A(z) = Ylk>m( Zk ~ Zm ) * s the Vandermonde determinant; the invariant 5-function 
can be represented as sum over irreducible representations 

s u{N) (n,v~ 1 ) = J2 d R x R (vn) (37) 

R 

where da is the dimension of the R'th representation and X R iP) is the character 

Xr(T exp i Adt) = Tr fl f exp (^i dt^ AjT^ , (3.8) 

where is a generator of U(N) in the R'th representation. 

It is well known [1] that in the case of free boundary conditions (i.e. ip(0) independent 
of <p((3)) the singlet partition function describes non-interacting fermions. Indeed after 
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integration over all angular variables only two Vandermonde determinants at the ends of 
the interval remain: A(z(0))A(z(f3)) and assure the antisymmetry of wave functions. In 
the case of the periodic boundary conditions, one should be more careful and use 5-function 
to match the eigenvalues: 

fW)) = [ n^(0)A 2 (^(0)) f(z(0))S(z(0) - z{(3)) ( 3. 9) 
J k 

Substituting eqs. (3.1), (3.6) and (3.8) in eq.(2.1) we obtain after integration over A(t) a 
functional integral over eigenvalues 

i r N 

' {v} fc=1 (3.10) 

where ^[ ^{^("^ ^ s a standard anti-symmetrizator, which appeared because of the 
skew-symmetry of the Vandermonde determinant: 

AiVz^V- 1 ) = (-l) v A(z(0)) (3.11) 

and Sq is the representation independent part of the action: 

S (z) = N f\t[^ + V(z)} (3.12) 
Jo z 

The T-ordering acts in the space of a representation R (r^rjl is the matrix product in this 
space). As it follows from eq. (3.2), the action does not contain a quadratic part for the 
diagonal elements of the field Aij. The integration over them gives rise to a selection rule 
for representations. It is shown in the Appendix A that it selects only such R for which 
the sum of all components of the highest weight, 



R = [mi, 777-2 • • • ™>n] 
m\ > m2 • • • > wtv, (3.13) 

equals to zero f: 

N 

]Tm fc = (3.14) 



fc=i 



f Let us note here that this class of allowed representations is much wider than the class of 
the self-conjugate ones (i.e. such that rriN-n = N — m n ) which was considered in the paper [3]. 
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We conclude from eq. (3.10) that Z N (fi, A) can be represented in the form of Gibbs 
partition function: 



Pr defined by eq. (A5) is a projector onto the subspace of all zero weight vectors in the 
space of representation R (i.e. the kernel space of the generators of the Cartan subalgebra). 
The sum in eq. (3.15) runs over all irreducible representations obeying the selection rule 
(3.14). Only for such representations this subspace is not empty. 

We see that unlike the simplest case of the trivial representation the eigenvalues cannot 
be considered here as non-interacting fermions (see Appendix A where transformational 
properties of ^/-functions under permutations of eigenvalues are described). 

A natural question arises whether, for higher representations, the same double scaling 
limit is applicable as for the trivial one, since the interaction among eigenvalues could 
completely change the critical properties of the system. However, we hope that, if the 
higher representations really correspond to the implementation of vortices, the critical 
properties of the world sheet of string are still the same as without vortices, at least 
far away from a vortex position. So, it is conceivable to adopt the same double scaling 
procedure as in [1] . It is based on the fact that the potential V(z) has a shape as shown in 
fig.l. It has a well filled by eigenvalues, which are kept there since the tunneling through 
the barrier is exponentially small with respect to its width. The decay of this system 
is not essential at least in every order of the (topological) expansion. The instability 
corresponding to the dominance of large graphs (thermodynamical limit) emerges when 
the eigenvalues near the Fermi level begin to spill over the top of the barrier, which can 
be achieved by tuning either A or N. 




(3.15) 



R 



where the Hamiltonians Hr are given by 




(3.16) 
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Fermi level 




Fig. 1: A typical potential V(z) for an eigenvalue z in eq.(3.12). The ground state 
is fermionic, i.e. first N quasi-stable levels in this unstable potential are occupied. The 
criticality appears when the Fermi level \x p touches the local maximum at z = zq. 

This shows that one has to concentrate oneself on the vicinity of the top which defines 
all the interesting physics. Let us choose new coordinates 

x i = VN(z i -z ) (3.17) 

where zq is the position of the top (zq = j for the cubic potential (2.2)). Then 

n^)) = -^ 2 +E Ck 



2N~ ' ^ N k / 2 

fc>3 



k x k 



(3.18) 



all the powers of X k 7 k > 3, are suppressed by the negative powers of N and can be omitted 
in the double scaling limit. Therefore, the essential part of the Hamiltonian (3.16) is 



N 



i a 2 i i 

2 x k\ + 2 ^ (r- - r V 



2 O.v'l 2 



(3.19) 



Of course, this Hamiltonian is ill defined because of the unstable quadratic potential 
and at some moment we have to remember about higher powers of x in eq. (3.18). The 
cubic term in eq. (3.18) provides an infinite wall at the distance 

(3.20) 



from the top of the potential. 
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As was noticed in ref . [1] , one can define a completely stable system having the same 
1/N-expansion by putting two symmetric walls at distances ±A from the top, i.e. our 
system can be stabilized by putting it into the N-dimensional hyper-cubic box {—A < 
Xi < A, i = 1, • • • N}. The convenience of the eigenvalue representation (3.19) is based on 
the fact that every eigenvalue has independent boundary condition. However, unlike the 
singlet representation, the eigenvalues are now interacting. 

On the other hand, it is obvious that eq. (3.19) can be viewed as the eigenvalue 
representation of the matrix upside-down oscillator with the Hamiltonian: 

%)=Y(^ + ^V) (3-21) 

where we have introduced the frequency ui by simple rescaling of the time. On the first 
glance, eq. (3.21) represents the system of N 2 non-interacting oscillators 

1 d 2 1 N d 2 

flM = -2S(5^)-2g(5^) ^ 

But the boundary conditions, which are determined now by higher terms of the type 

trcp 3 = y~](pij(pjk<Pkj (323) 

i,3,k 

mix all matrix elements and are not so simple as those for eigenvalues. 



4. Oscillatorial partition functions in different representations. 

The results of the previous section show that we can view one dimensional bosonic 
string theory as a theory of the upside down matrix oscillator, eq. (3.21), with U(N)- 
symmetric stabilizing walls (cut-off) at a large distance (~ A) from the top of the potential. 

Formally, if we forget for a while about the cut-off, this system is related to the 
standard oscillator with the stable matrix potential by simple analytical continuation to 
the imaginary frequency: 

uj -> iu (4.1) 

We might hope that physical quantities for the standard oscillator such as the partition 
function, Green's functions etc. already contain some information about the upside down 
oscillator. 
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But if we take, for example, the partition function for the standard oscillator, which 
is equal to the product of partition functions of all matrix elements: 

(1/2 \ N 2 

1_\ (4.2) 

where 

q = e-P« , (4.3) 

we see that it is too simple to be able to describe all the complexity of physics of 1- 
dimensional strings even after an analytical continuation. Apparently, the information 
about the cut-off is lost in (4.2), and we have to try more suitable quantities in order to 
get some information about the upside-down oscillator. 

A natural idea arises to classify all the states of the oscillator with respect to the 
irreducible representations of the U(N) rotations of the matrix coordinate: 

<p -> Q+ipQ (4.4) 

The corresponding wave functions *&r(<p) of these states should transform as matrix ele- 
ments in a given representation R. Being formulated as a function of eigenvalues of tp this 
wave function obeys the Schrodinger equation: 

H R y R (z) = EV R (z) (4.5) 

with H R given by eq.(3.19). 

The partition function in the R'th representation will be given by 

4") = Tre-f 3 "" (4.6) 

As we shall see in the next section, this quantity contains much more information 
about the continuation to the upside-down oscillator. The success of this approach is 
based on the fact that the cut-off wall for the upside-down oscillator is defined by higher 
U(N) invariant terms in the matrix potential (such as Xitrip 3 + \2trip 4 + ■ ■ •), and therefore 
it does not destroy the classification of states in representations: every wave in a given 
representation R being sent to the wall will be reflected with the same properties under 
U(N) transformations. Hence, we think that the analytical continuation procedure should 
be separate for every representation. 

Let us now work out the framework of finding Z R N ^ for the standard oscillator. 
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The most elegant way to do it is to define the so called twisted partition function: 

Z^(fl) = Tr(e-^e(Q)) = J dep a {E)e-^ E (4.7) 

where 6(0) is a twist operator which "rotates" the final states with respect to the initial 
ones by the U(N) matrix O and 



pn(E) = -ImJ2 [ d N \ 



N 2 v* E ,((p)y E ,(n+(pn) 



E-E' + ie 



(4.8) 



can be called the "twisted density of states" though it does not have even normal properties 
of a density (say, it is not positively defined). 
Then we can represent the eq. (4.7) as 



z( N \n) = J2d R z£ ) x R (n) (49) 

R 



and for the twisted density of states: 

Pq{E) = d R p R (E)x R (n) ( 410 ) 
R 

where the characters Xr(Q) obey the orthogonality relation: 

J (dQ)x Rl (n + )X R2 (n-uj) = 6 RlR2 X Rl (uj) (4.11) 



and hence 



and 



z R N) = J (dn)x R (n)z^ N \n) (4.12) 

p£\e) = J (dn) XR (n)Pn(E) (4.13) 

The simplest way to make these formulae more explicit is to use the Green function 
defined by the equation: 

&r tr( \w-^ )]G( ^' v ' )=0 (414) 

with the initial condition: 

G(0,<p,<f/) = 6l N *\<p-< f /) (4.15) 



15 



The solution is well known to be 



G(A tp, f') = /2 exp{ - |cth^ tr^ + ^) + i^MW)} (416) 

After the simultaneous change 

uj -> zcu (4.17) 

/? -> i/3 (4.18) 

eq.(4.16) can be viewed as the Green's function for the upside-down matrix oscillator. 

Eq.(4.16) appears to be a very useful starting point to many problems related to the 
one-dimensional strings. Let us, for example, demonstrate that the singlet sector of this 
theory, i.e. the sector bound to operators of the type 

JV 

O a (t) = tre a W = e QZfe(t) (4.19) 
k=i 

which are invariant functions of ip depending only on eigenvalues, can be described in terms 
of non-interacting fermions. 

Indeed, if we calculate the two-point function, D ai0l2 ((3), we obtain: 

D aiOC2 (j3) =« £re Ql ^°W Q2 <^ >= 

J d N \ J d N \'tre aiip tre a ^ 'G(/3, <p, <p') = 

2 (■ N 

(^)-e/ n^^xp{-|^? + ,?)} 

k,m i=l 

x e a ^+ a ^ x A 2 (z)A 2 (z') J (dQ)exp{ 7 ^J2 Q ab4^baZ a } 

a,b 

2 p N 

= £ / U"^m - \<*^i + *?) + £^.*u 

x A(z)A(z> aiZfe+a22 - (4.20) 

where we used the "angular" decomposition of the matrices <p = uj + zuj and 
ip' = {VLuj) + z' {VLuj) and the well-known formula [16] : 
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We see from eq. (4.20) that the initial and final states in terms of eigenvalues appear 
to be completely anti-symmetric due to the Vandermonde determinants A(z') and A(z). 
Hence, they can be represented as Slater determinants of the eigenfunctions of N oscillators 
(or upside-down oscillators): 

^nm2-n N (zi, ■•■z N ) = — L= det ^m(zj) (4.22) 

VM (*>j) 

So, we deal with free fermions. In principle, starting from here we can obtain all the 
results of the paper [11] for the two-point function in the double scaling regime. 
On the other hand, if we take the two point function of another type 

K&,=<%[O ai (0)O a M]> 

= I d N \d N \' t ^[e^e^']G{L3^) (4.23) 

we have a more complicated integral over U(N) group of the type 

(dn)Cl±Q ji exp{a tr(Q+z'Qz)} (4.24) 



Since the matrix element of the adjoint representation looks as 

= (0 + ),A, - ±515} (4.25) 

we expect that the Green's function (4.23) describes the propagation of states transforming 
under the adjoint representation of U(N). Indeed, one can show (see Appendix B) that 

K aia2 (P) = E(°l( eaiZm )( e "^ a *)r (e a2Zn ) |0> (4.26) 



m,n 



where 



A«)r="«E(-^-^)+ 



+ -S mn S 



r " j v§^-^) 2+ £(^-^) 2 ) fe-4) 2 (427) 



which coincides with eq. (3.16) if we take rf^ there. The simplest of these quantities is 

Ko(J3) = ^^-K QlQ2 (/3)| ai=Q2=0 =< tr[<p(0)<p(p)] > (4.28) 
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Fig. 2: A typical </? 3 Feynman graph for the adjoint propagator (4-28). It has a 
topology of a disc with two points on the boundary fixed. 

In the case of the matrix </? 3 -potential, Kq((3) will be defined by the diagrams of the 
type drawn in fig. 2. One can see that it represents an open string amplitude with free 
boundary conditions and two points on the boundary fixed in the t-space at and (5. 

On the other hand, this diagrammatic picture allows us to interpret the partition 
function (4.12) in the adjoint representation directly as describing the one- vortex-anti- 
vortex-pair sector of the model. 

Indeed, open string is the only intermediate state propagating in the adjoint represen- 
tation. For Z a dj we have the diagrammatic representation drawn in fig. 3. It looks similar 
to that of the correlator (4.28)) (fig. 2) but, instead of having "in" and "out" states at the 
ends, corresponding to some particular operators <p(0) and <p((3), it is continued period- 
ically looking now as being wrapped around a cylinder. In fact, there are two loops in 
the graph which wrap around the time circle of the length /?, in opposite directions with 
respect to each other (if we introduce, say, the clock- wise orientation for every loop). As 
we mentioned in section 2 it is exactly the definition of vortex and anti-vortex encircled 
by these two loops, correspondingly. All other loops are small in the t-space and do not 
contain vortices. 
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B 



Fig. 3: A typical Feynman graph for the adjoint partition function Z^ . It has 
the topology of a cylinder, which is wrapped around the circle of the radius (3/2tv. 

It is possible to generalize this picture to higher representations. 

From these considerations we draw the conclusion that higher representations have to 
be related in some way to open strings with the free (non-renormalized by some additional 
particles at the ends) boundary conditions and, hence, to bigger numbers of vortices. 

Now let us proceed with the partition functions. The untwisted partition function can 
be defined as: 

Z^(P) = J d"*<pG(/3,<p,<p) (4.29) 

Using eqs.(4.16) and (4.29) we can easily reproduce the trivial result (4.2). 
The twisted partition function is equal to 

Z( N \P,Q) = J d N2 cpG{p,cp,n + (pfl) (4.30) 

Here we simply have rotated the final states with respect to the initial ones by a matrix 

a 
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Let us first consider the trivial representation (singlet). Using eqs. (4.12) and (4.16) 
and diagonalizing again the matrix <p, we integrate over the angular variables by eq.(4.21) 
and obtain: 

r N 

fc=i 

x det ij [exp{-^-^z l z J }} (4.31) 

It is obvious from here that del in eq. (4.31) produces the expected fermionic statistics for 
{zk}- Say, the contribution of configurations with coinciding coordinates z m = z n for some 
m and n, is zero. Later we shall compute Z s i ng by other means and reproduce well known 
results for non-interacting fermions. The measure in eq.(4.30) is invariant under the group 
transformations ip — > uj + puj. We can use this invariance to show that Z( N \/3,fL) depends 
only on the eigenvalues of O. Indeed, we can use the ^-transformation of <p to diagonalize 
Q, since nothing depends on u: 

(uVtuj + ) kj = 5 kj e ie i (4.32) 

Let us perform now the gaussian integration over cp in eq.(4.30). We obtain the useful 
formula: f 



N 1 



k,m=l 

From eqs. (4.12), (4.30) and (4.33) we obtain: 

i r 27r N ria N 
Z «W = W II IT I A ^ I' *^ ie ) II h r a i -(a a J ( 4 - 34 ) 

where 

A(e*) = \{(e^-e^) (435) 

fc>m 

and the character Xr is given by the Weyl formula: 

X*(e'») = det ^g; ,m) (436) 



f We thank M. Douglas who showed us this formula. He claimed it was known long time ago 
but we did not find any references. 
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Using the complex variables 

6 = e ldk (4.37) 
eq. (4.34) can be written in the form of the contour integrals: 

y N 2 /2 - , N 



X 



N 

x J] (i-^-^r 1 (4-38) 

m,fc=l 

The integrals are taken over the unit circle, and all the poles inside this circle are concen- 
trated either at the origin or on the circle of the radius q. 

In principle, we can obtain by this formula using the direct computation some results 
for given R and N. However, we need the general formula for arbitrary N. 

Let us work out a more effective method for it. 

Using the "generating function for characters" [17]: 

N 1 

and the fact that the character is a homogeneous function: 

X R (pO=P^^ rnk X R (0 (4-40) 
where {m} are again the highest weight components, we can represent eq. (4.34) as: 

(4.41) 

where the integer numbers obeying the condition 

£ 1 >£ 2 --->e N (4.42) 
are related to the components {n^} of the highest weight of R' by the formula 

£ k = n k +N-k (4.43) 

and 

' J ° k=l 

21 



is the multiplicity of a representation R" in the tensor product R <8> R' of irreducible 
representations R and R' . For our purposes we need only the case R' = R" . It is important 
to note that only non-negative values of £k appear in eqs. (4.41), which makes the sum 
convergent. For the trivial (singlet) representation, we have Xsing(e* e ) = 1, and eqs. (4.41)- 
(4.44) give: 

oo oo oo 

4>)=E E ••• E ^-<* + »» = 

e N =oe N -i=i N +i £i=£2+i 

q N 2 /2 

= (l-q)(l- q 2)...(l- q N) ( 4 - 45 ) 

We have obtained the partition function of N fermions in the harmonic potential, as 
it was expected. The fermionic statistics appears here because of the ordering condition 
(4.42) for the components of the highest weight (see ref. [18] for some analogies). 

The formula (4.41) gives another check of the selection rule mentioned in the section 
3: Zr is nonzero only for representations R = {nik} with the sum of highest weight 
components equal to zero (eq.(3.14)). 

If the conditions 

JV 

ilk - n k+1 > y^(mj+ | m N |) = N | m N | (4.46) 
i=i 

(where {n^} are the components of the highest weight of R' in eq. (4.41)) are fulfilled for 
all 1 < k < N, then [17]: 

= 4 0) ( 4 - 4 7) 

where d^) is the dimension of the subspace of all zero weight vectors in the space of a 
representation R defined by eq. (A6) and in general 

M$aR> < d R (4-48) 
Eq. (4.47) means that, for (3 -> (or q = e"^ -> 1 ), 

d {0) 

Z R (3) ~ R (4.49) 

The general structure of the partition function (q) is such that 

4 w> (</) = -prw4>) c.50) 
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where Z^ g (q) is the singlet partition function defined by eq. (4.45) and P^\q) is a 
polynomial in q of a degree not higher than N ■ m, where m is the sum of positive (or 
negative) components of the highest weight of R. 
From eq. (4.49) we have: 

i#°(l) = 4 0) (4-51) 

Pft(q) can be regarded as the (/-analog of the multiplicity of the zero weight and is a 
known object in modern group theory f. They appear to be a particular case of the Kostka- 
Green-Foulkes polynomials [19]. We elaborated an effective recurrence algorithm suitable 
for us (see Appendix C for details) which enabled us to obtain for two infinite series of 
Young tables 

An = (n,0, ,-l,-l,...,-l) (4.52) 

and 

B n = (n,0, ,0,-n) (4.53) 

the following simple answers for the corresponding generating functions 

oo N-l 

V ( A N \x,q) = J2PA n \<l)x n = n( 1 + ^) ( 454 ) 



n=0 k=l 
N-l 

xq h 



n=0 fc=l q 

Eqs. (4.54) and (4.55) have the form of the grand canonical partition functions for fermions 
and bosons, respectively, in the N — 1 level system where x plays the role of the fugacity. 
Unfortunately, answers for other representations cannot be represented in a similar simple 
form but, nevertheless, one can obtain for representations having a small number of non- 
zero components of their highest weights formulae valid for arbitrary N. 



5. Partition functions for the upside-down matrix oscillator 

Now let us try to generalize the methods of the previous section to a more complicated 
case of the upside down oscillator. 

The major complification is the necessity to introduce a cut-off in order to stabilize 
this unstable system. Any cut-off decreases the high symmetry of the oscillator and makes 



f We are indebted to A.Kirillov for drawing our attention to this fact. 
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an exact solution impossible. The only tool in our hands is the analytical continuation in 
the frequency but we should be ready to confront with divergences and, hence, we need 
additional guesses in order to get physically reasonable answers. 

In order to classify the partition functions with respect to U(N) representations, we 
have to define again the twisted partition function. As it was noted in the section 4, we 
can formally obtain it changing 10 — > iu in the Green function G(f3, ip, ip'). Then, after the 
Gaussian integration over ip in eq. (4.30), we formally get the twisted partition function 
for the upside-down oscillator 

JV 

Z^(P,e) = e-^^ 11 [l-e-^+^-e™)]- 1 (5.1) 

k,m=l 

Instead of trying to give a meaning to eq.(5.1), let us make an attempt to obtain 
Z( N ^ (/3, 9) solving directly the wave equation 

HV((p) = eV((p) (5.2) 

for H taken from eq.(3.22). 

If we forget for a while about the boundary conditions (i.e. the cut-off) we can consider 
eq. (5.2) as a collection of independent wave equations for off-diagonal matrix elements: 

1 d 2 

~2 ( d<pi d(p*- + U VWW^V&ii) = tijVij&ij) ( 5 - 3 ) 
for i < j 

and diagonal matrix elements: 

1 d 2 

~2 K ~&p 2 +a;2 ^)*"(^") = e «*»(^«) ( 5 - 4 ) 
for z = l,2---iV 

If we introduce the parameterization (dropping for a while the indices (i,j)): 

^ = Vre ie (5.5) 

we get the following form of eq. (5.3): 

" 2 . r/^\2 e m 2 — 1-, n / r /?\ 
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where X is related to ^ as 

oo 

vjr(r,0) = v /f *rn(r)e ime (5.7) 

m= — oo 

Eq.(5.6) coincides with the Schrodinger equation for the hydrogen atom with the energy 
(f) 2 , the charge e/2 and the angular momentum t = l 7 "^" 1 . The imaginary values of 
cj in eqs. (5.3) and (5.4) correspond to the discrete spectrum of the hydrogen (negative 
energies), and we get: 

e=| u | (| m | +1 + 2&), fe = 0,l,2--- (5.8) 

The case of real u; in eq.(5.6) (the upside-down oscillator) corresponds to the continuous 
spectrum of the hydrogen, and we have to introduce some cut-off to define it. 

Again, as in the case of the infinite time interval [6] , we have strong reasons to believe 
that, for the exact definition of the spectrum, only the quasi-classical asymptotic of the 
wave function suffices: 

/ s 2 . / e , , n s 7r I m I — 1 I m I +1 ie A , _ 

X m (r) ~ - sm [r + -log(2r) - + ar,r(^^- + _) j (5.9) 

1 e 2 
for r > — (m 2 H -) 

2a; a; 2 

The Schrodinger equation (5.4) for the diagonal elements describes for an imaginary 
lo a collection of standard oscillators having the equidistant spectrum: 

e n =\u\(n+^), n = 0,l,2-.. (5.10) 

and for a real u the solution of (5.4) is a parabolic cylinder function having the quasi- 
classical asymptotics 

ty(x) ~ x~2 eX p| ± z[x 2 H logx+ 

+ argT(i^ + h} (5.12) 

2 

valid for a; 3> 

The problem now is that we are, in principle, not allowed to define the energy spectrum 
for every matrix element, cpij = y/fiJexpiOij and <pa = x, separately, since we have the 
U(N) invariant boundary condition, say: 

V(<p) = for tnp 2 > A 2 ~ iV (5.13) 
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which mixes all matrix elements. 

Strictly speaking, to solve this problem we have to consider first the whole wave 
function for all matrix elements, in the form 

{e}{m} i>j k 

and than find the coefficients C({e}) from the condition (5.13). But we are not able to 
find them directly. 

Instead of doing this we shall assume as in ref. [1] that the essential part of the 
spectrum is defined by the phase of reflection from the quadratic potential (last terms in 
the phases of the asymptotics (5.9) and (5.12) ) and not by the infinite wall, appearing in 
the conditions (5.13). 

To justify it let us note that every term in the expansion (5.11) has the quasi-classical 
asymptotics of the form 



exp ± % [tnp 2 + — ^2 € ij lo S + 



w .... . 
i>j 



+ argTi 1 "^^ + + £ ^gT^-f + t)] (5 .15) 

i>j " ~" k 



The first overall term in the phase of eq. (5.15), e ±ltrLp does not have any influence 
on the density of states, since it is invariant under U(N) transformations, ip — > fi+cpfi, as 
well as the boundary condition (5.13). 

The second logarithmic term, as we hope, will give only slow dependence on the cut- 
off A in the density of states, as it was in the simpler situation for the eigenvalues of the 
matrix <p in the singlet representation [4-7]. 

In other words, we expect that the density of spectrum for every matrix element 
can be defined from the phase quantization condition. Hence, from (5.9) we have for the 
off-diagonal degrees of freedom: 



A e , . 9 7r I m I — 1 _ / 1 m I +1 it x 



n = 0, 1,2- •• (5.16) 



which gives for the density of states p {e) = ^ 



„<~> (£)= __i_ Be ^ + ^I±l) + J_ i09A2 (5.17) 
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where 



^) = -J1(^-t-t 1 -t]- c 



(5.18) 



k=0 

is the ^-function of Euler; C is the Euler constant. 
For the diagonal degrees of freedom we obtain: 

A 2 + —logA 2 + argF(i— + ^) = irn 



u uj 2' 

n = 0, 1,2- •• (5.19) 

or, for the density of states Pd(e) = 

Pd(e) = - — Re^i- + I) + — log A 2 (5.20) 

TVUJ UJ 2 TVUJ 

Let us now calculate the twisted partition function for the upside-down matrix os- 
cillator. We shall start from the "twisted density of states" for the off-diagonal matrix 
elements. 

Since cpij = ^frije %Bii , the twist ip — > fi + </?fi results in 

$ij — > 9ij + 9i — 9j (5-21) 
and we conclude from eqs. (5.7) and (5.17), that 

oo 

p(9,e)= e* m V m) (e) 

m= — oo 

i °° imO( \m\+l , 7 \ 

= -E E 2 11! +<^)logA 2 (5.22) 

The cut-offs which we have introduced till now are not U(N) invariant. Therefore, we shall 
drop in what follows the last term with the periodic 5-function having in mind that the 
logarithmic divergency will appear on the next stages of the calculations and we will have 
to recover a cut-off again. 

The remaining double sum can be easily calculated (taking j = | m | +1 + 2k as a new 
variable) and the result is: 

sinh^(7r-#) 

p(0,e)= . V • a (5-23) 
sinh-7rsin6' 



It is remarkable that the twisted density of states does not contain a cut-off parameter 
for 9 7^ 27rm, m = 0, ±1, ±2 . . .. The twisted partition function of a non-diagonal matrix 
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element can be calculated taking both parts of the spectrum above and below the top of 
the potential, i.e. for positive as well as for negative e: 



+00 



Z (6,f3) = / dee- p *p(6,e) 



—00 



1 f +0 ° R sinhMix-6) 

For 

< e < 2tt (5.25) 
the integral in (5.24) is convergent. The calculation gives: 

Zo(e,P) = - a (5.26) 

cos [Ju — cos 

For 6 = and 2n we understand eq. (5.26) as an analytical continuation of eq. (5.24). 
If we rewrite eq. (5.26) as 

1 /2 P -i/3w 

(5.27) 



cos[3uj - cos# (1 - e-^+*))(l - e -Kfa>-0)} 

we find that this is quite similar to the factor 1 _ qei( i e k -e m ) in eq.(4.33): the differences is 
only in the change u — > iu. It means that we do not need to introduce a cut-off for the 
off-diagonal matrix elements. For the diagonal matrix elements ipa , the partition function 
does not depend on twist angles but the problem is that the integral over energies 

/oo poo 1 

de e-^p^e) = - de e'^-Re^i- + -) (5.28) 
-00 J-00 n u) 2 

is highly divergent and needs a cut-off procedure. However, we can formally deform the 
contour of integration from the real axis to the one shown in fig. 4 and expand ^-function 
in pole terms (eq.(5.18)). Then we can calculate the integral in eq. (5.28) term by term: 



h=0" w 2 + (^ + 2 ^ n =0 

1 

2sin/^/2 (5 - 29) 
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X 



X 5/2 LOJ 
X 3/2 tCJ 
X 1/2 tco 



-1/2 tCJ 
-3/2 
-5/2 tCt) 



Fzg. 4 -'The contour of integration for the diagonal matrix elements in the partition 
function (5.28). 

If we keep assuming that all N 2 degrees of freedom (pij for the upside down matrix 
oscillator are independent , we can write down the whole twisted partition function as a 
product of N 2 factors (5.26) and (5.29) 

Z {N) (0, P) = (2 sin <f) ~ N J] R 1 m r . (5.30) 

which coincides with eq.(5.1), and the partition function in a given representation R takes 
the form: 

1 r 2ir N m 

Z R N) ^ = ^l \{ d ^\^ e )\ 2 X R ^)Z^\e,{3) (5.31) 

It looks exactly the same as eq.(4.34) for the standard oscillator but with the change 
uj — > iu. This is an encouraging analogy for us if we remember that we used for the 
derivation of eq. (5.30) the true spectrum of the upside-down matrix oscillator. 

Let us also note that in eq. (5.31) there is no factorization to iV 2 independent degrees 
of freedom. Effectively, they interact and after integration over the angles only N effective 
degrees of freedom remain (e.g. iV fermions in the singlet sector). 
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Formaly, we can again write eq. (5.31) in the same form as eq. (4.41): 

Zg\p) = Y,M s ms e~^U^ (5.32) 

s 

with the multiplicity given by eq.(4.44). But one has now to give a meaning to the 

imaginary energy levels. 

For this purpose let us analyze the simplest case of the singlet representation for which 

M^ s = 1 (5.33) 



Hence, 



{e K >£ K+1 } 



(5.34) 



which can be reproduced from the result for the singlet representation (4.45) of the standard 
oscillator by change u> — > iu. 

How to give a meaning to eq.(5.34)? Let us pass from the micro canonical to the grand 
canonical ensemble. The singlet partition function takes the form 

oo 

oo 

= expJ2 l °9( 1 + e ^~ Mk+ ^ )] ) ( 5 - 35 ) 



k=0 



in the complete analogy with the well-known expression for the grand canonical partition 
function of fermions in the oscillatorial potential. As usually, the complex energy levels 
have to be understood as resonances [1,4-7]: 



/-t-oo i ^ k + - 

-oo 77 fc =0 ^u) + y K+ 2> 

/oo 
dep S ing(e) log(l + exp/?(/x - e)) (5.36) 
-oo 

where p S mg(e) coincides with eq.(5.20). 
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It is assumed that every integral in the sum is to be calculated as a pole contribution 
of the denominator in eq. (5.36) ignoring the singularities of log. In fact, at the lower limit, 
the integral in eq. (5.36) is logarithmically divergent and one has to cut the integration by 
some minimal energy eo(] eo |~ N), which is defined by the bottom of the original potential. 
Furthermore, the sum over resonances also diverges and needs to be regularized as well. 
Nevertheless, we have reproduced the correct result representing / sing as the partition 
function of N fermions in the upside-down oscillatorial potential [2] . Now, assuming that 
the general structure of answers is the same for all representations, we can proceed and 
get for the adjoint using eq. (CIO) the following formal expression: 



and, in the grand canonical form, 



oo 



Z^(3) = J2e^ N d adj Z^(P) = 



N=0 
2 



\^-\ -l\(A(f,,p)Z sing ( f ,-iu;,P) + B( f ,,p)) (5.38) 



where we have represented the dimension of the adjoint ,d a( y = N 2 — 1, as the differential 
operator and pulled it out of the sum over N. 



If we define that 



then 



„iu>(3 c — iujf3N 



B(/»,/J) = (5.41) 

In any case, this constant is determined by the first terms in the grand canonical expansion 
in eq. (5.38) and cannot influence universal properties of Z ac y 

Now using the same procedure which has led us from eq. (5.35) to eq. (5.36) we 
obtain: 



ld_ 



- ljA((i,/3)exp J dep ad} (e) log (l + e^" e )) (5.42) 
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where 

1 (z 3 1 

Padj(e) = Reifj(i- + -) + — logA (5.43) 

TVUJ UJ 2 TVUJ 

or 

1 1 

Padj(e) = p sing (e) - 2ivuj T (5.44) 

v a; / 4 

We see that the effective density of states in the adjoint representation differs from 
the singlet one by one missing Lorentzian (resonance). 

The factor A(fi, (3) could be rewritten with the same argumentation in the form 

/oo 
dep ad] (e) (5.45) 
-A 2 

and interpreted as the partition function of one particle decoupled from the others, even 
though we do not have any microscopic arguments supporting this interpretation. Never- 
theless, it seems to give just a no n- universal factor of the order exp(/3A 2 ) which does not 
affect universal behavior. Of course, this particle does interact with others at distances of 
the order A, but, being taken into account, this interaction produces only exponentially 
small in N terms and does not change the double scaling limit. 

Unfortunately, except for the adjoint representation, the procedure of the analytical 
continuation is ambiguous. If we collect together four representations of the next order 
(corresponding polynomials given by eqs.(Cll), (C12), (C13)) and implement the same 
procedure as led us to eq. (5.38), we shall get 

Z 2 (jjl,0) = J2 e^ N {2d A2 Z { /J(P) + d B2 Z^(P) + dcX N 2 \P)) 

N 



(l-e-^)(l-e- 2i ^)\ 2 \\(3dfiJ J\\f3dn 
1/1 9\Vl 9 A/1 3 



-Am) {m^Am-vs^-^ 
+ (r=^MW -'Mm) ~ 4|+ 

1/1 9\ 2 /l 9 \/l 9 



1//1 9\ 2 \ / / 1 d N 2 



-'AKm) - 4 (6 - 46) 
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where we have substituted d Aa = \(N 2 - 1){N 2 - 4); d B2 = jN 2 (N - 1)(N + 3); d C2 = 



Eq. (5.46) should contain, as we believe, contributions describing different configu- 
rations having the total vortex (anti-vortex) charge equal to 2. Moreover, there should 
be a contribution from two vortex pairs living on different surfaces (since we deal with 
the partition function rather than the free energy). In order to perform an analytical 
continuation, one should presumably separate all those terms and treat them differently. 

6. Physical results in the double scaling limit 

The results for the upside down oscillator given by eqs.(5.36) and (5.42) allow us to 
obtain the double scaling limit for the compactified one dimensional string (or for the 
XY-model coupled to gravity) in the sector with one vortex-anti-vortex pair. 

The partition function for the singlet representation eq. (5.36) was investigated in 
ref. [2]. It corresponds to the case of no vortices, when only nontrivial cycles of a two- 
dimensional manifold (given by a Feynman graph with non-trivial topology) can wrap 
around the time circle of length [3. 

Singlet partition function 

Let us repeat the calculation of / s ing(/?, fS) from ref. [2] for the consistency of content. 
We have to remember that the eq.(5.36) describes a filling by free fermions of both sides of 
the inverted oscillatorial potential. However, in the string perturbation theory which will 
correspond here to 1/ /(/-expansion we have to take into account only the states on one side 
of this potential, which corresponds to the perturbative vacuum of this string theory. This 
is reflected by an overall factor 1/2 in the following formula for the singlet free energy: 



where /J^ 3 (aO = A (\) + A 1 (\)fi + A 2 (A)/i 2 + • • • is the regular part of expansion of the 
grand canonical free energy in fx. Using an integral representation for p S m g (e) 



(we put uj = 1, fixing appropriately the scale of e) and integrating over e in the eq. (6.1), 
we obtain the following formula for the free enery of singlet states 



jN 2 (N + 1)(N — 3) and converted them into the differential operators. 




(6.1) 





(6.2) 
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Expanding / sing in and introducing a cut-off A to reguiarize the logarithmic divergences 
we get: 



\f log^/A 2 ) - 1(1 + (^) 2 ) log^/A 2 ) + fx^- 1 ^) 



2 r~ ovr~/ / 24 v ■ V p 

where we introduced polynomials in 4 



k=2 

(6.3) 



/.(« = (» -3)1 2-»£ (^r { ^" > -g.f" "^ff-' 1152 " 1 

^ /5 [2(fc - n)]![2nj! 

and i? m are the Bernoulli numbers. 

We do not have a direct interpretation of this grand canonical partition function as 
a sum over world sheets embedded into a time circle and classified by specific topologies. 
Such interpretation exists only for the canonical ensemble where we fix the variable N 
instead of \x. The relation between two ensembles is given by the integral transform: 

exp[/™ s ] = fdfieM-^N + r^^)] (6.4) 

where the integration contour encircles the point = 0. The 1/N expansion of f^ n9 (X) 
provides the string partition functions of given genera. 

In the double scaling limit which is the continuum limit of the string theory we have 
to extract only the lading singularity at A — > A c 0. To do it let us rewrite (6.4) in the form 

exp[f s ™ 9 } = j d/zexp^/^olog/io + f sin9 (fi)} (6.5) 
where we introduced the parameter /io by the formula 

^ologWA 2 ) = 2tt[Ai(A) - N] ~ NA 

with A = (A 2 - A 2 )/A 2 . 

The integral can be calculated in terms of the series in 1/fio which is equivalent to the 
1/N expansion. To generate this series we have to apply the saddle point approximation 



1 A c is defined in such a way that in the planar limit N — > oo when the temperature j3 1 is 
effectively zero (the spherical string amplitude does not feel the compactness of "time") we have 
A*(A C ) = 
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and then to take systematically into account the fluctuations around the saddle point. The 
saddle point for h(Ho) can be found from the equation: 

M„logWA 2 ) = ^^|^ (6-6) 
Solving it iteratively we obtain: 

(6.7) 



11 = HQ 



1 -i^te^ i - 1 H +0( i^ 5 j ) 



We see that in the double scaling limit when Ho is kept finite and A — > oo we can take 
H = Ho as a saddle point. 

Then we have to put h = h(Ho) + ^H and integrate over S/i in the Gaussian 
approximation. The contribution to f^ n9 from this Gaussian integration looks as 
— \ log log (//o/ A 2 ) + Q( iog( M „/A^) )- These terms should be neglected with the same ac- 
curacy. 

The net result of this calculation is that with the accuracy up to inverse logarithmic 
corrections we substitute the integral (6.5) by its saddle point value at \i = Ho an d get the 
following result for the canonical free energy of singlet states [2] : 



rsvng 



p [^logWA 2 ) - ^(i + (^) 2 ) ioswa 2 ) +jrH - 2ik - 1) fM 

k=2 



24 v v (3 



(6.8) 



Note the change of sign of the first term (spherical free energy) with respect to the similar 
term in (6.3). 

As it was noticed in ref. [2] eqs. (6. 2), (6. 3) and (6.8) possess a T-duality symmetry 
with respect to the inversion of the radius of compactification: 

(3 2n (3 

This important property remains to be true in every order of the topological (1/N) expan- 
sion (we recall that h ~ A" | lo ^ A in the planar limit). 



Adjoint partition function 
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The analytical continuation of the partition function for the adjoint representation 
(eq. (5.42)) looks ambiguous but nevertheless we assume that its universal part can be 
written as follows 



1 d 2 

Zad o = (^p" 1 ) ex P 0*> P) ( 6 - 9 ) 

where fadjin, ft) is determined similarly to f s in g (fJ>, P) but with the density of states given 
by eq. (5.45) 

1 1/2 

Padj(e) = Asing(e) - - g2 - ^ 4 (6.10) 

Hence, we have 

5 = / adj - /sing = ^ J de(p ad] (e) - p sing (e)) log (1 + e^" £ )) (6.11) 
Using the representation 

/•OO 

: / dxe tX£ ~ zx (6.12) 



e 2 + z 2 



Re 



and integrating in eq. (6.11) over e we obtain for the third derivative of 5 the following 
representation 

OS 1 f°° , sin(/ix) sinh| 



a„ 2 ; dx 55£f (613) 



And the topological expansion of 5 takes the form 



log(/l/A2) __L (1 _i ( | rt 



V u- 2m f-l) m f2m - IV 2" 2m V f 27r l 2n ( 22n - 2 )^2n 
2^ llj l** 1M [2(m-n)]!!(2n) 



m=2 n=0 

From eq. (6.14) we conclude that 



(6.14) 



/adj-Zsing^^log^ (6.15) 

Since \x ~ A • N in the first approximation (and A 2 ~ N), we reproduce in eq. 
(6.15) the result (1.1) from the paper [3]: the spectrum of angular (vortex) excitations is 
separated from the singlet (vortex-free) spectrum by a logarithmically big gap. Note that 
the expression (6.14) is not self-T-dual anymore. 

This fact can be qualitatively explained from the comparison of the sum-over-surfaces 
picture for the correlator eq. (4.20) propagating the singlet states and the correlator (4.28) 
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propagating the adjoint states. The first one is given by the sum over all Feynman graphs 
with two vertices fixed (on a graph and in the target space). These two vertices are 
arbitrary and can belong to any two different loops or the same loop on a graph. 

In the case of the correlator (4.28) one has to pick up two vertices belonging to the 
same loop on a graph, which is only a part of (positively defined) contributions of the 
correlator (4.20). Hence 

°~ D aia2 ((3) \ ai=a2=0 = D (/3) > K (/3) (6.16) 



If we take for the definition of the gap the large (3 asymptotics: 

D (/3) ~ exp( 

^sing 0) (6-17) 

ifo(,3)~exp(-m adj -/3) (6.18) 
we come from eq. (6.16) to the conclusion that 

TO ad j > m sing (6.19) 

Hence, it is not surprising that 

^sing- ,. 1 A , ^ (6.20) 

I log A I A^O v ' 

and 

WadjH log A | -> +oo (6.21) 

As we demonstrated in the section 4, higher representations correspond to the multi- 
loop amplitudes of open strings with free boundary conditions on the edges of the world 
sheet. We have no parameter (like the mass of a particle at the end-points of the open 
string) to adjust in order renormalize the "boundary tension" and to make the character- 
istic size of the boundary macroscopic. Eqs. (6.20) and (6.21) show that the characteristic 
boundaries have very short lengths, and the contribution of few vortices is negligible in 
comparison with the contributions of vortex-free configurations. 

However, the vortices have a considerable entropy, which causes at some (3 the 
Berezinski-Kosterlitz-Thouless phase transition [13,14]. Knowing the partition function 
for one vortex-anti- vortex pair, we can calculate the critical value, (3kt, as follows. From 
the eqs. (3.15), (6.8) and (6.15) we obtain in the spherical approximation in the canonical 
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(fixed N ensemble) the following partition function of the vortex anti-vortex pair in the 
c = 1 string theory [3] : 



Z N (P, A) ~ Z sing (1 + N 2 e f °*- f ' in ° +■■■) 

~ exp ( ^iV 2 ^A 2 [log(A) + const A^" 2 + •••]) (6.22) 
\2ir 2 J 

We see that the second term (as well as the contributions with more vortices and anti- 
vortices) is irrelevant with respect to the first one, describing the usual scaling of the c = 1 
string with no vortices, when (3 > 4ir. However, it becomes more important than the 
first term for (3 < 4n and the dilute gas approximation for vortex-antivotex pairs used in 
the eq. (6.22) is not valid any more. The world sheet of the string will be immediately 
densely populated by the strongly interacting vortex plasma. It is a typical picture of the 
Berezinski-Kosterlitz-Touless phase transition. 

This phase transition occurs in our model at the same inverse temperature 

Pkt = 4tt (6.23) 

as in the model of planar rotors (XY model) on the plane. Indeed, the answer (6.22) 
should be compared with the configuration integral of the vortex-anti-vortex pair on the 
plane [13,14] (in the same normalization): 

Q= [d 2 x [ d 2 yexp-^log |x ~/ 12 ~V 2 ~£ (6.24) 
J J 2tt r z Q 

where V is the volume of the space. 



7. A possible approach to D+l-dimensional bosonic string 

A natural generalization of the matrix model describing one-dimensional strings to the 
case of D+l-dimensional bosonic strings is the following D+l-dimensional scalar matrix 
field theory: 

Z N (P,\) = J V N2 <p(x 7 t)exp-Ntr J alt J d D x[^(d t "<pd fl <p) + V(<p)] (7.1) 

with the potential V((p) given by eq.(4.3) or of a more general form. 
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The analogy with the string theory is again based on the Feynman graphs: we obtain 
for the free energy = -j^logZffi the following diagrammatic expansion 

oo oo „ 

f N ((3) = N- 29 E A " E / dD+lx i ■ ■ ■ d ° +lx - II D ( x * - x i) (7.2) 

where 

D ^ = W^S dD ^ e ' iri, 7^ < 7 - 3 > 

As long as we have the superconvergibility of the integrals in eq. (7.2) (it is true, say, 
for D < D c = 5 for </? 3 graphs) we may hope that the particular microscopic definition of 
D{x) is irrelevant, and for big enough k's in eq. (7.2) we have defined the D + l-dimensional 
bosonic string field theory. 

Unfortunately, the theory (7.1) is much more complicated than the theory (2.1). We 
cannot, for example, obtain a representation similar to the Hamiltonians (3.16) for the 
eigenvalues of the matrix field: the corresponding "connection" = fi+c^fi (compare 
with eq.(3.3)), obeys now the more complicated constraint, d^A v — d v A^ + [A^, A v \ = 0, 
and the integration over A^ is highly nontrivial. The importance of the angular variables 
for C = D + 1 > 1 physics is reflected in the fact that the eigenvalues may be now not the 
only and even not the most important degrees of freedom. 

On the other hand, if we still hope that the physics at C > 1 bears some universal 
features, we may admit that the particular shape of the potential V(tp) should not be 
important. All the results for C < 1 show that only singular points like local maximum or 
inflection points are essential for the continuum limit of the corresponding string theory, 
and all the rest in the shape of the potential after an appropriate rescaling serves as a 
cut-off. 

This brings to mind the idea that it may be sufficient to consider the most generic 
unstable potential: 

V(<p) = ~<f? (7-4) 

Forgetting for a while about the rest of the potential, we can use the solvable Gaussian 
theory 

Z N = J V N2 (p(x)exp - ytr J d D+1 x[{d^) 2 +m 2 ^ 2 } (7.5) 

in order to get the results for the unstable theory (7.4) by means of the analytical contin- 
uation 

m — > im (7-6) 
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and an appropriate introduction of a cut-off. If we will be unable to find this procedure 
so the whole approach based on the matrix model (7.1) does not describe any universal 
bosonic string theory for c > 1. 

Again as in the case c = 1, the Gaussian theory (7.5) is too poor to be a starting 
point for an analytical continuation, and we have to consider the states obeying some 
particular representation of the U(N) group (since the potential trV(<p) is still invariant 
under <p — > Q + </?Q transformation and so is the cut-off procedure). 

It is natural to hope that the ground state of the system is again the singlet under 
the group transformations; therefore, we have to extract only the singlet part of eq. (7.5) 
and then to perform the analytical continuation. 

In order to find the partition function for the singlet states we introduce again the 
twisted partition function: we compactify the time dimension on the circle of a length (3 
(t G (0, (3) again) and for the remaining D directions we use the momentum representation: 

Z%(0,Q) = J V N2 <p p (t) expytr ^ dt J d D p[ip 2 - (p 2 + m 2 )^ 2 ] 

<p p (o) = n+<p p (/3)n (7.7) 

where we have introduced the U(N) twisted periodic initial conditions. 

For fixed p every matrix element can be considered as an independent oscillator with 
the effective frequency 

oo p = vV + m 2 (7.8) 
Hence, in the complete analogy with eqs. (4.33) or (5.30) we can write: 

N 

Zg\/3 t d) = exp- E nWk-Oj) (7.9) 
k>j=i 

where 

F(P, 9) = L D J d D p log [cosh^ vV + m 2 ~ cos9] (7.8) 

L is an infrared cut-off (the size of the box where the system is put in the D dimensional 
target space), and the twist matrix O is chosen in the diagonal form 

= diag(e i9l ,e i9a ,---e i9N ) (7.9) 
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Then the recipe to find the ground state E of our model is to calculate the singlet 
partition function: 

r 2ir r 2ir N m 

4>) = / • • • / n P- 1 ^ ie ) i 2 z n d) ^ o) (T.io) 

Jo Jo k=1 ^ 

to perform the analytical continuation (7.6), introducing the cut-off on the way, which 
should defined the dependence on the cosmological constant, and to find the asymptotics 

Of course, all these steps do not seem to be trivial. One of many obstacles to overcome 
is the highly divergent integral in the eq. (7.8). one has to learn how to deal with these 
divergences. But what is hopeful here is the fact that the integral in eq. (7.10) goes over 
only N (and not iV 2 ) variables and is in principle of a saddle point type in the large N 
limit. 

Another way to proceed is to try to use eq. (4.39) in order to obtain a representation 
similar to eq. (4.41): 

z { R ] m = n 1 A{ei9) |2 XR(ei9) n n e F w»-w 

J ° k=l 71 {p}k<j 

= x) e "^^= i( ^ ,p+i) /wxii(^n^p(^)xs P (^) 

S p " { P } 

Sp 

where £k,p are defined analogously to eq. (4.43) and attached to each point of the p space. 

^ s sum over all representations at all points. In the continuum limit it becomes 
a fermionic path integral. M^^J j plays the role analogous to the corresponding factor 
in eq. (4.41) making the system highly non-trivial. It gives rise to the same selection rule 
(3.14) for allowed representations R as in the one-dimensional case. 

We hope to come back to this problem in the future and in this paper we propose 
eqs. (7.9), (7.10) as a definition of a D+l-dimensional bosonic string theory directly in the 
continuous limit, as it was discussed in the introduction. 

Another interesting question is the role of D-dimensional partition functions in higher 
representations. They might describe vortex-like excitations (i.e. monopole) which might 
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be quite important for the physics of strings. It is not clear whether they are separated by 
a gap from the singlet vacuum. 

8. Discussion. 

The aim of this paper was to demonstrate the possibility to extract physical infor- 
mation about ID bosonic string by means of the analytical continuation of quantities 
known for the standard matrix oscillator. It was shown that not only the singlet (with 
respect to the U(N) symmetry of the model) but also the adjoint partition functions can 
be obtained in the double scaling limit with the help of this procedure. It was shown 
that the adjoint partition function describes the one-vortex-anti-vortex-pair sector of the 
XF-model coupled to gravity. It enabled us to calculate the critical temperature for the 
Berezinski-Kosterlitz-Thouless phase transition. 

However, the method of analytical continuation obviously suffers from ambiguities 
which appear because the high symmetry of the matrix oscillator cannot be preserved 
in the upside-down case. Any cut-off imposed in order to stabilize the system decreases 
this symmetry. Therefore we had to guess a dependence on the cut-off using close anal- 
ogy between the singlet and adjoint partition functions. However, in the case of higher 
representations, one apparently needs some microscopic information in order to make the 
continuation unambiguous. To begin with, one needs e.g. the decomposition of the sum 
over all representations with respect to the contributions of particular vortex-anti-vortex 
combinations. It may be done matching the partition functions with propagators of non- 
singlet states as it was done in this paper for the adjoint representation. Another useful 
trick would be to introduce a regularization similar to one proposed in ref. [11] in the 
singlet case. 

It is tempting to interpret the non-local matrix variable 



as an operator creating a vortex of the unit charge. It would be interesting to include the 
angular degrees of freedom in the collective field representation [22] as well as eigenvalues. 

Although, knowing the partition function for one vortex-anti-vortex pair, one can 
correctly estimate the position of the Berezinski-Kosterlitz-Thouless phase transition, in 
order to get a solution at the critical point, one should sum over all or, maybe, over the 




(8.1) 
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most important representations. This problem is far from being solved. However, it seems 
that properties of the XF-model interacting with 2D quantum gravity are quite similar to 
those of the XY-mode\ on the plane: for f3 > 4tt one has the dilute-gas-of-dipoles phase, 
where a vortex and an anti-vortex are tied together forming a dipole; and for (3 < 4-7T one 
has the plasma phase, where interactions between vortices are screened and an emerging 
finite correlation length leads to the change of the matter central charge from c = 1 to 
c = 0. 

The model under investigation is the first example of a matrix model where the angular 
degrees of freedom (and not just the eigenvalues) of the matrix field play an important role. 
Now, it is clear that the angular degrees of freedom describe vortices and hence they are 
of a great concern in many physical applications. They become crucial when the central 
charge is bigger than 1. Apparently, one has to pay much more attention to them in order 
to understand the nature of this phase of the string theory from the matrix model point 
of view. 
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Appendix A. 

In this appendix we derive the selection rule (3.14) for representation contributing to 
the sum (3.15) and discuss the behavior of eigenf unctions of the Hamiltonian (3.16) with 
respect to permutations of the coordinates zi. The Fourier decomposition of ^/-functions of 
the original matrix Schroedinger equation gives rise to the expansion of the Green function 
(defined in the case of the harmonical oscillator by eqs. (4.14) and (4.15)) 



where A and fi are eigenvalues of NxN matrices X and Y ( X = fiiAO+,y = 02/^0^) 



means that eq.(Al) is invariant under left and right shifts of U by a diagonal unitary 
matrices: Uu — > Ukie l( ~ 0k ~ (fl \ By definition, 




(A.l) 



{R} a, 6=1 





(A.3) 



and, from the invariance of the Haar measure under this symmetry, 



J{dU)G{\U»U+)D%{U) \[ d ^D a Ae ie ) JJ \[ d -^D Vh {e-^) 



(A.4) 



P a a'9a'b'(\ ^)Pb'b 



where Pji 



; 27r n 



k=l 2tt 



D a b(e lB ) is a projector , since 



ab 



(A.5) 



The dimension of the subspace it projects on to is equal to 




(A.6) 
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Under permutations of A(A — > pXp 1 ) the Fourier coefficients behave as follows 

g5,ipXp-\ fi) = J (dUMU+pXp-^^DKiU) 
= D*,(p)g5 b (\ tf i) 

and analogously for p 

g^ppp- 1 ) = J {dU)G{\Upw- 1 U+)D% t {U) = 
gn,(\,p)DK b (p- 1 ) = D b R bl (p)g^(X, f i) 



(A.7) 



(A.8) 



It is clear from eqs. (A7) and (A8) that g^ b (X,p,) is a tensor operator with respect to the 
permutations, i.e. the Weyl group. It means that properties of the Green function (or the 
^/-functions) with respect to the permutations are determined by the matrix elements of 
the Weyl group, which are non-trivial in general. For self-consistency the subspace defined 
by eqs.(A4) - (A6) has to be invariant under all permutations of A and p. It is not difficult 
to show that it is the subspace of all zero-weight vectors in the space of a representation 
which form an orbit with respect to the Weyl group. 

All vectors in the space of a representation can be numbered by the Gel'fand-Zetlin 
patterns: 

m 1 Nm 2 N rn N -iNm N N 



mi,N-i m2,N-i mN-2,N-i rriN-i,N-i 



(A.9) 



mi3 m 2 3 to 33 
mn 

The rows in eq. (A9) are the highest weights in the following sequence of embeddings [19] 

U(N)DU(N-1)D DU(2)DU(1) (A.10) 

where U(fc) is spanned by the first k rows and columns of the U(N) matrix in the fundamen- 
tal representation. The first row in eq. (A9) coincides with the U(N) highest weight and is 
fixed. Other numbers are restricted by the single condition : mk,i > rrik,i-i > rrik+ij- The 
character of U(l) is equal to exp(imn^i) and, hence, in order to get a non-zero answer in 
eq. (A6) one has to take mn = 0. From the identity U(/c)=U(l)xSU(/c), it follows that 
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X u {k){ e%tp ) = X S u(k)( elip ) ex P« Ylj=i m j,k{ ( P\ + ( P2 + ■ ■ ■ + l Pk) and one has to take the trivial 
representation for all such U(l) subgroups: 

k 

rfc = E mn > fc = (A.ll) 
i=i 

Eq (All) has to be fulfilled for all k < N and corresponding vectors will have the zero 
weight (by definition, the weight of a vector is equal to [ri, — r±, . . . , tn — rjv-i]). When 
k = N ', eq. (All) gives the selection rule (3.14). 



Appendix B. 

Using the decomposition (3.1) and the eqs.(3.2), (3.3), we can represent (4.23) (in the 
functional integral formalism (2.1), (2. 2), but with free boundary conditions for ip(0) and 
<P(P)) as 

N 

K ai aM= E (e ai ^ (0) |(O(0)O+(/3)) m jV 22 »^) (B.l) 

m,n=l 

Note that, for the free boundary conditions, the constraint (3.5) is absent due to 
the linear measure for the hermitian matrix field, A(t), with the Gaussian weight (see 
eq.(3.2)). This functional integral is ultra-local in time, and one has only to take care of 
time ordering. The Vandermonde determinants of the z-variables cancel one another after 
the integration over A(t) such that only two of them remain at the end points (see [1]) 

As the result we get 

N.N poo 2 2 

K ai aM= E /n^W e "" tr/ - dt(4Zfe+7Vl/(Zfe)) ] A W--))A(,(oo)) 

m,n=l k = l 2) 
x e a l2m (0)+a 22 „(/3) Ap exp j__L f P dtQ(z)} 

V 2A" J Q 

where the matrix Q(z) is of the form 

~ / I 1 \ 5 m 5 n 

Q( z )™ = 5 mn 5 l3 V + V ^ - 2 . 1 \ 2 (B.3) 

where the double indices (i,m) and (j,n) represent ones in the space of the adjoint repre- 
sentations. The unitarity condition, 00+ = i", leads to the equations 

N N 

E®? = E^? = °> ( R4 ) 

i=l j=l 
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hence the number of independent indices is equal to the dimension of the adjoint, = 
N 2 -l. 

In the Hamiltonian language eqs. (B2), (B3) are equivalent to eqs. (4.26), (4.27) 

2 

(where it was already taken V(z) = —\)- 
Appendix C. 

One has to calculate the sum over ordered integer numbers of the type 

S fl (g)= 2^ M ms q^i=i ni ^ CA j 

ni>n i + 1 

where rii are the components of the highest weight of a representation S. The coefficients 
Mft^g can be found by the following graphical rules [21]. Let us draw the vertical line and 
place the Young tableau of a representation R in such a way that the number of boxes to 
the right of the line equals to the empty space to the left. After that one has to add all 
boxes of the tableau of a representation S to the first tableau (R) so that to reproduce the 
second one (S) to the right of the line according to the following algorithm. One has to 
place all boxes from the first row of the second tableau in different columns of the first one 
so that they have formed a possible Young tableau. Then one has to proceed with boxes 
from the second row but with the restriction that above and to the right of each box in an 
obtained tableau the number of them must not exceed the number of boxes from the first 
row. Then one has to repeat this procedure for all rows. It is convenient to attach labels 
to all boxes of the second tableau, boxes from the same row having the same labels. Then 
MjL, s will be equal to the number of all possible combinations of labels. The structure of 
labeling is not fixed only in columns in which there are boxes of the first Young tableau. 
Therefore, one can consider every labeling in these columns as a separate sum over rii-s 
contributing to eq.(Cl) which takes in this case the form 

3*(«) = ££9 Ei = in< (C2) 

{L} {nt} 

where 1S ^ ne sum over au possible labelings; J2{ ni } * s tne sum over N integer numbers 

corresponding to all possible tableaux compatible with a particular labeling from the first 
sum. All numbers run to the infinity but from below the sum is restricted by a condition 
more complicated than the simple ordering. Nevertheless this sum gives just 

q a(L) E sing (q) (C.3) 
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with an integer a(L), hence eq.(4.50) is true and 



P R (q) = J2 ( i a{L) 

{L} 



(C.4) 



a(L) = J2 kA k 



(C.5) 



k=i 



where = mining — rik+i), (un+i = 0) are minimal differences between the indices in 
the second sum in eq.(C2). 

If all labels less than k have been already placed in some way, one has a configuration 
which is equivalent to the case when the number of rows is equal to N — k and placing 
k-labels one gets a number of configurations with the number of rows less by 1. It gives 
rise to recurrence relations between contributions to Pii(q) from different tables with the 
decreasing number of rows. For example, for the sequence of representations defined by 
eq. (4.52) there are only two ways to place fc-labels: 



k k 



k k 



\k\k\ 



k k 



a k - a k+l + <? a k+l 

and PA n (q) = a^K Eq. (4.54) is a direct consequence of eq. (C6). For representations 
(4.53), the recurrence relation is a little bit more complicated: 
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&l n) = £«*&£Ti (C7) 

and PB n (q) = (q). It is clear that eq. (C6) describes fermions (no more than one 
particle can fill the /c'th level) and eq. (C7) describes bosons (the number of particles 
are not restricted). One can proceeds in the same way for other representations but the 
structure of answers will be the more complicated the more non-trivial Young tableaux 
are. Nevertheless, they can be always interpreted as statistical sums for n particles in the 
equidistant level system, where n is the number of positive (or negative) components of a 
highest weight. For example, for the sequence of representations 

,0,-l,...,-l) (C.8) 

N v ' 

n 

one gets the recurrence relation: 



C n = (l,...,l,0, 
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i=l 

which describes fermions but with "wrong" weights: the weight of a fermion on the k-th 
level is equal to q k ~i where j = if the (k — l)-th level is empty and j = k — 1 otherwise. 
The form of the polynomials for several simplest representations are: i) for the adjoint 

N 

ii) for 4 representation with n = 2 § 

q 3 (l-q N - 2 )(l-q N - 1 ) 



Pa 



(l-q)(l-q 2 

= T-> T\ l u - ii J 



Pt 



C-2 



q 2 {l-q N - 3 ){l-q N ) 



(l-g)(l- ? 2) 
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